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3 ABSTRACT 

We present the first results from the X-Shooter Lens Survey (XLENS): an analysis of 
the massive early-type galaxy SDSS J1148+1930 at redshift z = 0.444. We combine 
its extended kinematic profile - derived from spectra obtained with X-Shooter on ESO 
Very Large Telescope - with strong gravitational lensing and multi-color information 
derived from Sloan Digital Sky Survey (SDSS) images. Our main results are (i) the 
luminosity- weighted stellar velocity dispersion is (<x*)0$ R e &) — 352 ± 10 ± 16kms _1 , 
extracted from a rectangular aperture of 1'.'8 x l'/6 centered on the galaxy, more 
accurate and considerably lower than a previously published value of ~ 450 km s -1 ; 
Q_i (ii) a single-component (stellar plus dark) mass model of the lens galaxy yields a 

q logarithmic total-density slope of 7' = 1.72lp;Qg (68 per cent CL; p tot oc r~ 7 ) within 

5— 1 a projected radius of ~ 2'.' 16; (hi) the projected stellar mass fraction, derived solely 

CZ) from the lensing and dynamical data, is /*(< Re) = O.IE^qqq (68 per cent CL) 

inside the Einstein radius for a Hernquist stellar profile and no anisotropy. The dark- 
matter fraction inside the effective radius /dm(< R c b) — 0.60+QQg ± 0.1 (68 per 
cent CL), where the latter error is systematic; (iv) based on the SDSS colors, we 
find /*,saip(< Re) = 0.17 ± 0.06 for a Salpeter Initial Mass Function (IMF) and 

CO /*.Chab(< Re) = 0.07±0.02 for a Chabrier IMF. The lensing and dynamics constraints 

t" ■ on the stellar mass fraction agree well with those independently derived from the SDSS 

I" colors for a Salpeter IMF, which is preferred over a Chabrier IMF at variance with 

standard results for lower mass galaxies. Dwarf-rich IMFs in the lower mass range of 

(T) 0.1-0.7M Q , with a ^ 3 (with dN/dM oc M~ a ) - such as that recently suggested for 

massive early- type galaxies with a = 3 in the mass range 0.1 — 1M & - are excluded at 

y—( the > 90 per cent CL. and in some cases violate the total lensing-derived mass limit. 

^ — 1 We conclude that this very massive early-type galaxy is dark-matter dominated inside 

one effective radius, consistent with the trend recently found from massive SLACS 

• i-h (Sloan Lens ACS) galaxies, with a total density slope shallower than isothermal and 

a IMF normalization consistent with Salpeter. 

H 
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1 INTRODUCTION 

Understanding the evolution and the internal structure of 
massive early-type galaxies (ETG), as well as their stel- 
lar and dark matter distributions, is crucial in order to 
fully comprehend the processes in hierarchical galaxy for- 
mation (e.g. White & Rees 1978; Davis et al. 1985; Frenk 
et al. 1985). In this context, the relationship between bary- 
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onic matter, which dominates astrophysical observables, and 
dark matter (DM), which dominates most of the dynamics 
during galaxy formation is particularly important. 

To unravel some of the these issues and paint a more 
robust physical picture, enormous effort has been afforded 
to study the relative contributions of baryonic, dark matter 
and black hole constituents of ETGs through stellar dynam- 
ical tracers, X-ray studies, and gravitational lensing (e.g. 
Fabbiano 1989; Mould et al. 1990; Saglia, Bertin & Sti- 
avelli 1992; Bertin et al. 1994; Franx, van Gorkom & de 
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Zeeuw 1994; Carollo et al. 1995; Arnaboldi et al. 1996; Rix 
et al. 1997; Matsushita et al. 1998; Loewenstein & White 
1999; Gerhard et al. 2001; Seljak 2002; Borriello, Salucci & 
Danese 2003; Romanowsky et al. 2003; Treu & Koopmans 
2004; Cappellari & Emsellem 2004; Cappellari et al. 2006; 
Thomas et al. 2007; Czoske et al. 2008: Czoske, Barnabe & 
2008; Auger et al. 2010; Treu 2010). 

The picture emerging over the last decades from studies 
of the inner regions of early- type galaxies, where baryonic 
and dark matter are both present, is that to first order the 
total mass density profile inside few effective radii can be 
well-described by a power-law form close to an isothermal 
profile with 7 ~ 2 for p tot — r -7 , although with a ~ 10 per 
cent intrinsic scatter in the density profile (e.g. Gerhardt et 
al. 2001; Treu & Koopmans 2004; Koopmans et al. 2006, 
2009; Auger et al. 2010, Barnabe et al. 2008, 2011). The 
dark-matter density profile in the same region, however, is 
far less well constrained although seems consistent with a 
density slope 7dm ~ 1.3 (e.g. Treu & Koopmans 2004). 

In addition, while the innermost regions of early-type 
galaxies are expected to be dominated by the stellar mass 
component, the dark-matter mass component is usually 
found to play a non-negligible role, with mass fractions rang- 
ing from 10 to 40 per cent of the total mass within one 
effective radius (e.g. Gerhard et al. 2001; Cappellari et al. 
2006; Gavazzi et al. 2007; Weijmans et al. 2008). Even more 
recently, observations, as well as theoretical studies based 
on stellar population and dynamical models (e.g. Bullock 
et al. 2001, Padmanabhan et al. 2004), indicate that, for a 
constant IMF, the dark matter fraction in the internal re- 
gion increases monotonically with the mass of the galaxy 
(e.g. Zaritsky et al. 2006; Auger et al. 2010), a trend that is 
more conspicuous in the case of slow-rotator ellipticals (Tor- 
tora et al. 2009). On the other hand, it has also been shown 
that the luminous stellar mass-to-light ratio scales with the 
luminous mass of the system (Grillo et al 2010). However, 
information from stellar kinematics is limited to the central 
regions (a few effective radii) because of the lack of bright 
kinematic tracers at large radii, and the total mass deter- 
mination suffers from the well-known degeneracy between 
the mass density profile of the galaxy and the anisotropy 
of its stellar velocity dispersion tensor (Binney & Mamon 
1982). Higher-order velocity moments, which potentially al- 
low one to disentangle this degeneracy by providing addi- 
tional constraints (Gerhard 1993; van der Marel & Franx 
1993), can only be measured with sufficient accuracy in the 
inner parts of nearby galaxies with current instruments. De- 
spite great progress, it still remains difficult to separate the 
stellar and dark-matter components, mostly due to a still 
relatively poor understanding of the precise shape of the 
stellar IMF and its associated stellar mass-to-light ratio. Un- 
certainties related to the latter can easily lead to a factor of 
a few uncertainty in the inferred stellar mass. 

Understanding the stellar IMF in massive early-type 
galaxies is a key open issue with a broad range of astro- 
physical implications. Although it is commonly assumed 
that the IMF is universal and independent of cosmic time, 
several authors have suggested that the IMF might indeed 
evolve (Dave 2008; van Dokkum 2008) or depend on the 
the stellar mass of the system (e.g. Worthey et al. 1992; 
Trager et al. 2000; Treu et al. 2010; Graves & Faber 2010, 
Auger et al. 2010). Recently, van Dokkum & Conroy (2010) 




Figure 1. SDSS color composite image of the sky region around 
the Cosmic Horseshoe (the zoom is 16" X 16"). 



have suggested that low-mass stars (< 0.3 M©) could be far 
more dominant in massive early-type galaxies than previ- 
ously thought. This could imply that much of the increase 
in the mass-to-light ratio of galaxies with galaxy mass is due 
to a changing stellar IMF rather than an increasing dark- 
matter fraction, consistent with the suggestions of Treu et 
al. 2010 and Auger et al. 2010. 

Valuable additional information on distant early-type 
galaxies, besides their kinematics, can be provided by gravi- 
tational lensing. Strong gravitational lensing provides a very 
useful tool to investigate the total mass content of the lens 
(Maoz & Rix 1993; Rusin & Ma 2001; Rusin & Kochanek 
2005; Dye & Warren 2005; Brewer & Lewis 2006b, 2008) and 
to place constraints on the inner mass distribution of lens 
galaxies at redshift out to z ~ 1. Unfortunately, the mass- 
sheet and the related mass-profile (Wucknitz 2002) degen- 
eracies do not always allow one to accurately determine the 
slope of the galaxy density profile and above all to unam- 
biguously disentangle the luminous and dark matter com- 
ponents. 

Gravitational lensing and stellar dynamics are partic- 
ularly effective to break many of these degeneracies, when 
they are applied in combination in the analysis of the in- 
ternal structure of distant early-type galaxies. The two ap- 
proaches are complementary and allow, when combined self- 
consistently, to robustly determine the mass profile of the 
lens galaxy (Koopmans & Treu 2002, 2003; Treu & Koop- 
mans 2002, 2003, 2004; Jiang & Kochanek 2007; Barnabe 
& Koopmans 2007; Czoske et al. 2008; Czoske, Barnabe & 
2008; Koopmans et al. 2009; Grillo et al. 2010; Treu 2010; 
Barnabe et al. 2011). Disentangling stellar and dark matter 
remains difficult, but the uncertainties have been reduced 
from factors of a few to far less than this (Treu et al. 2010). 

The X-Shooter Lens Survey (XLENS) aims to take the 
next step and spectroscopically observe a sample of lens 
galaxies from the Sloan Lens ACS Survey (SLACS, Bolton 
et al. 2006) with oetg 2 s 250 km s -1 using the X-Shooter in- 
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Table 1. Properties of the cosmic Horseshoe 1 . 





Parameter 


Values 


Lens 


RA 


llh 48m 33.15s 


Galaxy 


Dec 


19° 30' 03'.'5 




Redshift 


0.444 




Effective radii 


1'.'96 ± 0'.'02 




9L 


(20.84 ± 0.06) mag 




TL 


(19.00 ± 0.02) mag 




ih 


(18.22 ±0.01) mag 




zl 


(17.75 ±0.04) mag 




Axis ratio, g 


0.8 ±0.1 


Source 


Redshift 2 


2.38115 ± 0.00012 




Star formation rate 


~ lOOMQyr" 1 




Dynamical mass 


M vir ~ 10 10 Mq 


Ring 


Diameter 


10'.' 2 




Length 


~ 300° 




UL 


21.6 mag 




9l 


20.1 mag 



lL 19.7 mag 

Mass enclosed 3 (5.02 ± 0.09) X 10 12 Mq 

1 Belokurov et al. (2007) measured the redshift of the source 
to be z = 2.379. We find a systematic shift that brings the 
source redshift to be z = 2.3811 in agreement with Quidcr 
et al. (2009). 

2 The mass within the Einstein radius or, more precisely, 
within the ring diameter, is taken from Dye et al. (2008). 

s Parameters obtained from images taken with the 2.5 m 
Isaac Newton Telescope (INT). Magnitudes are taken from 
SDSS DR7. See Belokurov et al. (2007) 

strument . With this data we intend to further disentangle 
the stellar and dark-matter content of the galaxies, through 
combined lensing, dynamical and spectroscopic stellar popu- 
lation studies. With multi-band HST images in hand, we are 
able to obtain more precise dark-matter mass fractions than 
ever before, in order to ultimately correlate these with ETG 
mass and compare with theoretical galaxy formation mod- 
els. By combining the lensing and dynamical results with 
stellar population models we will plan to constrain the nor- 
malization and shape of the stellar initial mass function. 

In this paper we present the result of a pilot program 
of the XLENS project: A study of the "Comic Horseshoe" 
(SDSSJ 1148±1930), an almost complete Einstein ring with a 
diameter of ~ If/.' 2 around a very massive early-type galaxy 
at z = 0.444. The system was discovered by Belokurov et 
al. (2007) in the Sloan Digital Sky Survey Data Release 5 
(DR5, Adelman-McCarthy et al., 2007) . The source is a 
star-forming galaxy at z = 2.381 (Dye et al. 2008; Quider 
et al. 2009). Properties and characteristics of the Cosmic 
Horseshoe are listed in Table 1 and a composite SDSS image 
of the system is shown in Figure [l] 

The paper is organized as follows: In Section 2, we 
present the observations and data reduction. In Section 3, we 
discuss our kinematic analysis. In Section 4, we discuss the 
luminous and dark matter distributions of the lens galaxy. 

1 X-Shooter is a powerful broad-band (3000-25000 A), three- 
armed medium-resolution spectrograph on the VLT (D'Odorico 
et al. 2006); 22 h GTO time on three systems (084.A-0289 and 
087.A-0620) and 40^ GO time were awarded (P086.A-0312) to 
the XLENS project (PI: Koopmans) 



We summarize our findings and we present our conclusions 
in Section 5. We assume Hq = 70kms _1 Mpc -1 , fi m = 0.3 
and Q,a = 0.7 throughout the paper. 



2 OBSERVATIONS AND DATA REDUCTION 

X-Shooter observations of SDSSJ1148+1930 were carried 
out during a GTO run between 17-24 March 2010 in slit 
mode 2 , splitting the beam over three arms: UVB (R=3300 
with 1'.'6 slit); VIS (R=5400, with l'.'5 slit); and NIR 
(R=3300 with 1'.'5 slit), covering a wavelength range from 
3000 to 25000 A simultaneously. The 11" long slit was cen- 
tered on the galaxy with a position angle (PA) of 163°. 
The latter minimizes contamination from the source and 
leaves enough sky region to facilitate accurate sky subtrac- 
tion. Two Observation Blocks (OBs) were not used because 
of bad seeing and/or an incorrect positioning of the slit on 
the object. The total exposure time on target for each arm 
is ~ 7389 sec and the typical seeing is ~ 0'.'6. Standard 
calibration frames were obtained during daytime after the 
corresponding OB. A summary of the observing blocks is 
given in Table 2. 

Data reduction was done using the ESO X-Shooter 
pipeline vl.2.1 (Goldoni et al. 2006) and the Gasgano data 
file organiser developed by ESO. The pipeline reduction uses 
calibration spectra, taken during the commissioning run, for 
bias subtraction and flat-fielding of the raw spectra. Cosmic 
rays are removed using LACosmic (van Dokkum 2001). For 
each arm, the orders are extracted and rectified in wave- 
length space using a wavelength solution previously obtained 
from the calibration frames. The resulting rectified orders 
are shifted and coadded to obtain the final two-dimensional 
(2D) spectrum. We extract a one dimensional spectrum (ID) 
from the resulting 2D merged spectrum, using an IDL code 
that uses the optimal-extraction algorithm of Home (1986). 
It also produces the corresponding error file and bad pixel 
map. The final signal-to-noise ratio in the UVB±VIS spec- 
trum is ~ 7 per pixel. No telluric correction was applied, so 
that prominent atmospheric absorption bands can still be 
seen in the final spectrum (Fig[2|. 

Because the near-infrared spectrum suffers seriously 
from sky-line residuals using the current pipeline, we limit 
our analysis to the UVB- VIS region of the spectrum and 
defer a full analysis of the infrared data as well as our spec- 
troscopic stellar-population analysis to future work. In this 
paper, we base our stellar mass determinations solely on the 
broad-band colors from the SDSS plus the stellar kinematic 
and lensing data. 



3 STELLAR KINEMATICS 

We measure the luminosity-weighted velocity dispersion 
(LOSVD) of the lens galaxy from the final ID UVB-VIS 
spectrum using the Penalized Pixel Fitting (PPxF) code of 
Cappellari & Emsellem (2004). PPxF determines the com- 
bination of stellar templates which, when convolved with 
an appropriate line-of-sight-velocity distribution, best re- 
produces the galaxy spectrum. The best-fitting parameters 

2 P084.A-0289(A); PI: Koopmans 
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Figure 2. The final luminosity-weighted UVB-VIS ID X-Shooter spectrum extracted from a rectangular aperture of 1"8 X 1"6 centered 
on the galaxy (see text). Spectral features the lens galaxy are marked. Telluric absorption lines have not been removed from the spectrum 
nor has the spectrum been flux calibrated. The boxes indicate parts of spectrum that are affected by sky or telluric lines. 



Table 2. Observational details of the Cosmic Horseshoe 



Observation 
Block 



Date 



Exp. time 
(sec) 



Seeing 



200200337 


3-17-2010 


3 


X 


821(UVB 


- VIS) 


0'.'62 






3 


x 


3 x 274(NIR) 




200200343 


3-17-2010 


3 


x 


821(UVB 


- VIS) 


0'.'56 






3 


X 


3 x 274(NIR) 




200200436 


3-19-2010 


3 


X 


821(UVB 


- VIS) 


Of/66 






3 


X 


3 x 274(NIR) 





of the LOSVD are determined by minimizing a x penalty 
function. The best fit provides the mean velocity and the 
velocity dispersion (v and a, respectively), plus their uncer- 
tainties. The S/N ratio of the data is inadequate to measure 
the higher order Gauss-Hermite moments ft 3 and h^, which 
quantify the asymmetric and symmetric departures of the 
LOSVD from a pure Gaussian (related to the skewness and 
kurtosis respectively). PPxF also allows the user to mask 
noisy or bad regions of the galaxy spectrum. We also perform 
an iterative sigma-clipping in order to clean the spectrum of 
residual bad pixels, sky lines and cosmic rays. We focus on 
absorption lines between 3500-5500 A (including Ca K and 
H, G4300, H/3, Mgt, and some Fe lines). To minimize errors 
due to mismatch between the resolution of the templates and 
the galaxy spectrum, we use X-Shooter stellar spectra ob- 
tained as part of the X-Shooter Stellar Library (XSL) survey 
(Trager et al. 2011 in prep.), with similar instrumental reso- 
lution (for the galaxy spectrum we use a l'.'5 slit, correspond- 
ing to (Jinstr ~ 25kms _1 , while for the stellar templates the 
slit width is 0'.'7, corresponding to o"i nstr ~ 12kms _1 ). An 
excellent fit is obtained with a weighted linear combination 
of a Kl giant template (57 per cent) and a G2 star template 
(43 per cent). The selected region of our galaxy spectrum 



used for the fit and the corresponding best fit stellar tem- 
plate are shown in Figure [3] 



3.1 Luminosity- Weighted Kinematics 

The measured luminosity-weighted velocity dispersion for 
the central aperture of 1" 8x1'.' 6 is (a) = 352±10±16kms _1 . 
The formal error on the dispersion includes both the random 
error contribution and the systematic uncertainties due to 
spectral range differences, template mismatch and contin- 
uum fitting as discussed below. Figure [3] shows the best 
fitting template of the PPxF routine superimposed on the 
galaxy spectrum as well as the residuals of the fit. 

As a first test of the accuracy of our measurements, we 
use the more heterogeneous MILES 4 stellar template library 
(Sanchez-Blazquez, et al., 2006). We select 100 stars (F, G, 
K, M) in the range 3525-7500 A, with 2.3 A FWHM spectral 
resolution. The measured luminosity-weighted stellar veloc- 
ity dispersion of (a) = 358 ± 31kms , after instrumental 
correction, is consistent with the above estimate, based on 
XSL templates, but has larger errors due to the lower resolu- 
tion of the MILES library. As a second test, we fit templates 
to two different spectral regions. We find slightly different 
results between the blue and the red part of the spectrum, 
but always consistent within 2a (not including systematics) . 



3 The only velocity dispersion value previously published for the 
Cosmic Horseshoe is from Belokurov et al. (2007). They perform 
fits of Gaussian line profiles to the Ca H and K absorption lines 
from a much lower resolution spectrum (FWHM~ 10 A). They 
find a higher velocity dispersion estimation of 430 i 50kms _1 , 
inconsistent with our high-quality data. 

4 http:/ /www. iac.es/proyecto/miles/pages/stellar- 
libraries / miles-library.php 
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Figure 3. Top panel: Logarithmically rebinned UVB+VIS galaxy spectrum (black) and logarithmically rebinned best-fit template 
(green), both in restframe wavelength. Bottom panel: Residuals from the fit. Bad-pixels excluded from the fitting procedure (sky line, 
telluric lines) are shown in red. See the text for more information. 



The scatter in these fits is used to estimate additional sys- 
tematic uncertainties related to template mismatches and 
spectra coverage. 

3.2 Spatially- Resolved Kinematics 

To preserve the spatially-resolved kinematic information, we 
define seven spatially- varying apertures (with adequate S/N 
ratio) along the radial direction and we sum the signal within 
each aperture. Apertures are defined to be larger than the 
seeing, in order to have independent kinematics measure- 
ments for each aperture. The stellar rotation velocity and ve- 
locity dispersion are measured in each aperture using PPxF, 
as described above. Again, we use different spectral regions, 
excluding the most prominent telluric lines in the VIS range 
and a range of seven XSL stellar templates (G, K and M 
stars). The uncertainties on the inferred kinematics are esti- 
mated by adding in quadrature the formal uncertainty given 
by PPxF and the scatter in the results for different templates 
and spectral regions. Details of the aperture sizes and the 
kinematic profiles are listed in Table 3. The rotation and the 
velocity dispersion profiles are shown is Figure [4] 

We find an almost flat velocity dispersion profile be- 



yond the effective radius. The weighted average value of 
344±25 kms -1 is consistent within the formal error with the 
luminosity weighted value for an aperture of 1.8" x 1.6" (see 
Table 3), as expected. The velocity profile shows some mild 
rotation, with a projected rotation velocity of ~ 140 km s -1 
at one effective radius. Because the effective dispersion 
(» r 2 ms = y/v 2 + a 2 ; see Cappellari 2008) is well within the 
errors on a, the effect of rotation can be neglected in our 
spherical Jeans analysis and we will ignore rotation in the 
remainder of this paper. 



4 STELLAR AND DARK-MATTER 

Here we derive the slope of the total density profile and cal- 
culate the fraction of dark matter inside one effective radius. 



4.1 The Galaxy Mass Model 

To derive the stellar mass inside the Einstein radius, we fol- 
low Treu & Koopmans (2002), Koopmans & Treu (2003), 
Treu & Koopmans (2004). We model the mass distribution 
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Figure 4. Kinematic profiles of the Cosmic Horseshoe. Rotation (bottom) and velocity dispersion (top). Grey line shows the light- 
weighted value for the central aperture. 



as a superposition of the stellar and the dark matter compo- 
nents. For the stellar mass distribution we use two different 
spherical models, described by the equation: 



Table 3. Spatially resolved kinematics of the Cosmic Horseshoe. 



p,(r) = 



(3 - n)M*r« 
4 7tT n (r + r*) 4 " 



(1) 



The Hernquist (1990) luminosity-density profile has n — 1, 
and the Jaffe (1983) luminosity-density profile has n — 2. 
The total stellar mass is M* and the break-radius is r*. 

We model the dark matter distribution as a generalized 
spherical NFW profile, using a break radius of r& = 50kpc, 
is typical of massive ETGs (although the result is insensitive 
to its exact value) and an outer slope of 7 out = 3 (Gavazzi et 
al. 2007; Auger et al. 2010b), following the prediction from 
numerical simulations of dark matter halos (Ghigna et al., 
2000): 

Pd,0 



pDM { 



(r/ri,)-T[l + (r/n) 2 ]CT- 3 )/ 2 



(2) 



The break radius and the density scale {pd,o) determine the 
virial mass of the dark matter halo (Bullock et al., 2001). 

For completeness, we also use a single mass compo- 
nent model where the stellar and dark-matter mass den- 
sities add to a power-law with an effective density slope 7' 
for ptot oc r 1 and the stars are treated as trace particles 
(see Koopmans & Treu 2003). 

The most accurately known constraint from the lens 
mass model is the mass inside the Einstein radius. We use 
M E = (5.02 ± 0.09) x 10 12 M© inside E = 5"l derived by 
Dye et al. (2008). The mass components must add to Me 
within the Einstein radius. The error in radius is folded into 
the error in mass so that a fixed radius can be used (the 
mass and Einstein radius are coupled in the modeling). We 
use the average of the two effective radii from the g and 
i band images reported in Belokurov et al. (2007), R e s = 
1.96 ± 0.02 arcsec (Table 1). The effective radius is uniquely 
related to r* (see Hernquist 1990 and Jaffe 1983). 

We also assume a constant orbital anisotropy parameter 
ft and allow it to range between 0.0 and 0.5, which is typical 
for massive early-type galaxies (e.g. Gerhard et al. 2001). 
The parameters of the model without a strong prior are then 



Aperture 


Aperture 


V 


a 


center (") 


dimension(") 


(kms- 1 ) 


(kms" 1 ) 


0.00, 0.00 


1.80 x 1.60 


0+15 


352 ± 10 


-2.16,0.00 


0.80 x 1.60 


-100 ± 100 


350 ± 100 


-1.36,0.00 


0.80 x 1.60 


-80 ± 100 


311 ± 76 


-0.64,0.00 


0.64 x 1.60 


-9 ± 25 


341 ± 26 


0.00,0.00 


0.64 x 1.60 


0+12 


332 ± 16 


+0.64,0.00 


0.64 x 1.60 


62 ± 18 


360 ± 25 


+1.36,0.00 


0.80 x 1.60 


77 ± 80 


350 ± 100 


+2.16,0.00 


0.80 x 1.60 


180 ± 100 


410 ± 100 



the stellar mass (M*) and the dark-matter density slope (7), 
or only 7' for the single-component model. 

We subsequently solve the spherical Jean equations and 
compare the models to the kinematic data, taking the aper- 
ture sizes and seeing into account. We find an effective den- 
sity slope of 



7 



, 7o +0.05 
l - ' z -0.06 



(68percent C.L.) 



when marginalizing over f3 = 0.0-0.5, shallower than isother- 
mal (see Koopmans et al. 2006 and 2009). This slightly low 
value for the logarithmic total-density slope may suggest 
that this object can be a group, or a small cluster of galax- 
ies, where the overall efficiency of converting gas into stars 
is lower, and where typically the overall mass density pro- 
file in the corresponding region is shallower than isothermal 
(e.g. Newman et al. 2011). Changing the luminosity-density 
profile also does not change the final logarithmic slope sig- 
nificantly. The more interesting case of the two component 
model will now be discussed. 



4.2 The Stellar Mass Fraction from Lensing and 
Kinematic Constraints 

To derive the stellar mass fraction inside the Einstein ra- 
dius, we create a densely sampled grid of likelihood values 
by comparing the kinematic profiles to the data for a pro- 
jected stellar mass fraction [/*(< Re)] within the Einstein 
radius ranging between and 1 and a dark-matter density 
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slope (7) ranging between 0.0 and 2.0. We assume flat pri- 
ors on both quantities and marginalize over 7 to derive the 
probability density function of /*(< Re)- The results for 
the Hernquist and Jaffe profiles are shown in Figure 5 for 
/3 — and a more extreme case of radial anisotropy with 
P = 0.5. Taking as reference the best fit parameters of the 
Sersic profile computed in Dye et al.(2008): n = 5.40 ±0.04, 
r = 3.9 ±0.1, Li/2 = 61.2 ±0.4), with the form: 

L = L 1/2 exp{~B{n)[{r/T f- n - 1]} (3) 

we conclude that Hernquist and Jaffe profiles fit well the 
observed luminosity profile of the galaxy within the effective 
radius (with a slight preference for the Hernquist model). 
We find a fraction of stellar mass within the Einstein radius 
for the two luminosity profiles of: 

ifiq = 0.19to:SS and /j* F = 0.13lg;jg for f3 = 

and 

/Sq = 0.13±g;8r and /j* F = O.lll^s for = 0.5 

We note that in projection, the lens galaxy is already fully 
dark-matter dominated inside ~ 2.5 effective radii. For com- 
parison with previous work, we also derive the dark-matter 
fraction inside the effective radius and find 

/ D M(<-Reff) =0.60l^±0.1, 

[68 per cent confidence level (CL)] for /3 = 0, including a sys- 
tematic uncertainty of 0.1. The random error is based on the 
marginalized probability densities shown in Figure 5 and a 
systematic error is included based on the maximum range of 
possible dark-matter density slopes. Although the latter is 
constrained by the models, we extrapolate inward from _Re 
to J? e ff, where this slope could be somewhat different. We 
note though that the systematic error is rather conservative. 
The dark-matter fraction increases by ~ 0.1 for /3 = 0.5 and 
the difference between the Hernquist and Jaffe profiles is 
negligible (by construction, since they both contain equal 
fractions of mass inside that radius). This high dark- matter 
fraction inside the effective radius is consistent with the re- 
sult found in Auger et al. (2010) for SLACS systems and is 
consistent with the dark matter fraction within the effective 
radii beeing a monotonically-increasing function of galaxy 
mass. 

4.3 Stellar Mass Fraction from Stellar Population 
Constraints 

We independently calculate the projected stellar mass frac- 
tion inside the Einstein radius using stellar population syn- 
thesis models and the SDSS colors (see Table 1) of the lens 
galaxy. A comparison between this fraction and that de- 
rived from lensing and stellar kinematics provides a direct 
constraint on the stellar IMF (see also e.g. Treu et al. 2010; 
Auger et al. 2010). 

We use GALAXEV, a library of evolutionary stellar 
population synthesis models computed using the isochrone 
synthesis code of Bruzual & Chariot (2003). This code allows 
one to compute the spectral evolution of stellar populations 
for a wide range of ages and metallicities. Here we use the 
2007 version of GALEXEV kindly provided by S. Chariot 
(a version commonly known as 'CB07'). We use six Simple 



Stellar Populations (SSP) models computed with a Salpeter 
IMF or a Chabrier IMF for a range of metallicities from 
Z = 0.0001 to Z — 0.05, all of them with lower mass cutoff 
of Mi ow =0.1 Mq and upper mass cutoff of M up = 100 Mq. 
The spectral resolution is 3 A across the wavelength range 
of 3200-9500 A. 

We compute the spectral evolution of the stellar popu- 
lation and the redshift dependence of colors and magnitudes 
in the SDSS filters g, r, i, and z for all SSP models and for 
a range of ages and star formation histories (SFH). For the 
SFH model, we use an exponentially declining star formation 
rate (SFR) with time-scales r = 0.1 to 0.4 Gyr and as ex- 
treme cases an instantaneous burst with r = or a constant 
SFR. We compute the redshift evolution of the galaxy col- 
ors and magnitudes for a start of star formation between 12 
to 5 Gyr, corresponding to a formation redshift of between 
roughly z ~ 5 and ~ 0.5. By comparing the model magni- 
tudes to the SDSS magnitudes, we subsequently determine 
a grid (as function of age, metallicity and SFH duration) of 
likelihood values for each model as well as the total stellar 
mass. We do this for both the Salpeter and Chabrier IMFs. 

We use the standard Bayesian approach, as outlined 
in Auger et al. (2009), to determine the posterior proba- 
bility distribution function and the marginalized errors on 
the total stellar mass of the galaxy, assuming flat priors on 
all parameters in logarithmic space (e.g. a prior oc 1/r for 
t). The latter assumption is not critical, but given that the 
time-scale of the SFH, etc., are unknown a priori, this prior 
is a better description of our ignorance before making the 
observations (and modeling). From the resulting cumulative 
probability function, we calculate the median of the total 
stellar mass and its 68 per cent confidence interval for both 
IMFs. Using the observed brightness profile of the galaxy 
modeled as a Hernquist, Jaffe or deVaucoulers profile (the 
precise choice is not critical), we determine the fraction of 
the stellar mass (assuming a constant M/L ratio) within the 
Einstein radius. The results of this analysis are listed in Ta- 
ble]^ where we report the inferred total stellar mass fraction 
for each IMF as well as that from lensing and kinematics. 

We find that, for the range of assumed luminous pro- 
files and anisotropies, the stellar mass fraction from lens- 
ing agrees remarkably well with that from stellar population 
modeling assuming a Salpeter IMF (Fig. [5). This result fur- 
ther supports the results found by Treu et al. (2010) that 
the IMF of massive early-type lens galaxies is well-matched 
by an evolved Salpeter IMF. They found a tentative trend of 
the IMF slope with galaxy velocity dispersion, with a 'light' 
IMF such as a Chabrier IMF is inappropriate for systems 
with a ^ 200km s -1 . This trend may imply a non-universal 
IMF, perhaps dependent on metallicity, age, or abundance 
ratios of the stellar populations. Alternatively, it may imply 
non-universal dark matter halos with inner density slope 
increasing with velocity dispersion. While the degeneracy 
between the two interpretations cannot be broken without 
additional information, Treu et al. data imply that mas- 
sive early-type galaxies cannot have both a universal IMF 
and universal dark matter halos. This is confirmed by the 
expanded analysis of Auger et al. 2010, who added weak- 
lensing data to constrain the halo model, finding that the 
data for massive galaxies are inconsistent with a Chabrier 
universal IMF for a range of realistic halo profiles, including 
various recipes to account for baryonic effects. 
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Figure 5. Marginalised PDF for the Hernquist model (red lines) and the Jaffe model (blue lines) from our lensing and dynamics model. 
Solid lines are for an anisotropy parameter of = 0.0, and dashed lines for fi = 0.5. Black arrows show the internal stellar mass fractions 
with different IMFs from our stellar population analysis. Black lines represent PDFs of Van Dokkum & Conroy IMF and of the most 
extreme bottom-heavy IMFs. In the labels, the first number represents the break radius, the second is the IMF slope at lower masses 
(m ^ Tnbreak) and the last number is the slope at higher masses (fix to Salpeter slope). We exclude the Chabrier IMF and the a = —3 
IMFs at a confidence level of ^ 90 percent and the extreme bottom-heavy (a = —3.5) IMFs at a confidence level of 99 percent, while 
the Salpeter IMF is fully consistent with lensing data. 



5 BOTTOM-HEAVY IMFS 

Recently, van Dokkum & Conroy (2010) have suggested the 
presence of a large number of stars with masses ^ 0.3Mq in 
the central regions of early-type galaxies. By measuring the 
strength of the Nal doublet and the Wing-Ford molecular 
FeH band in the spectra of eight of the most luminous and 
massives early-type galaxies in the Virgo and Coma clus- 
ters, they infer that the IMF for these systems might even 
be steeper than Salpeter, with a slope as steep as a = —3. 
They also test even more 'bottom-heavy' with a = —3.5 and 
'bottom-light' (dwarf deficient) IMF but they find the best 
fit between stellar population synthesis models and spec- 
trum around the Nal doublet with the a — —3 IMFs, al- 
though the uncertainties are large and Salpeter cannot be 
excluded with high-confidence.. 

Here, we investigate this claim by assuming a broken 
power-law IMF, with the Salpeter slope in the high-mass 
regime that dominates the SDSS magnitides (i.e. changes in 
the IMF below this break do not affect the stellar-population 



analysis carried out above) and a steeper profile in the low- 
mass range. We test three different values of the break point 
in the mass function: mt re ak = 0.3,0.5 and 0.7Mq, respec- 
tively. Changing the IMF for stars with M sC mbreak does 
not change the SDSS colors because > 90 per cent of the 
light in the spectral region we studied here is coming from 
stars with M ^ 0.7 Mq. On the other hand, it dramatically 
changes the total stellar mass of the system, because stars 
with masses of 0.1 — O.3M0 can contribute at least 60 per 
cent of the stellar mass for these bottom-heavy IMFs. 

We calculate the change in total stellar mass arising 
from the change in the slope of the IMF: 



AM = 



JO.ll 



dN 




dN 




dm 


IMF 


dm 


Salp 



m dm. 



(4) 



where 
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Table 4. The stellar mass fractions within the Einstein Radius derived from lensing plus stellar kinematics and from stellar population 
synthesis models, respectively. The likelihood ratio compares maximum a-posteriori difference between the SSP model and the no- 
difference (null-hypothesis) model with the isotropic Hernquist model, which has the highest inferred stellar mass fraction. The IMF 
slope is assumed to be a = —2.35 beyond the break rribrcak m the mass function. The probabilities in between parentheses represent for 
each IMF the fraction probability that a given PDF (Pa) matches the max probability. 



Lensing and Kinematic Model 


Anisotropy 


Stellar Mass Fraction 




Hernquist 


p = 0.0 
P = 0.5 


1Q+ - 04 
u - iJ -0.09 

1S+ - 04 
u - ±o -0.07 




Jaffe 


p = 0.0 
P = 0.5 


n n +0.03 
u - ic> -0.05 

11+ - 02 
u - iJ --0.05 




Stellar Population Model 


^■brcak 

(Me) 


Stellar Mass Fraction 


21n(P max /P ) 
(HQ versus SSP) 


Chabrier 
Salpeter 

van Dokkum & Conroy a = —3.0 
van Dokkum & Conroy a = —3.0 
van Dokkum & Conroy a = —3.0 
Bottom- heavy a = —3.5 
Bottom- heavy a = —3.5 
Bottom- heavy a = —3.5 


0.3 
0.5 
0.7 
0.3 
0.5 
0.7 


0.07 ± 0.02 
0.17 ± 0.06 
0.30± 0.11 
0.39 ±0.15 
0.51 ± 0.18 
0.48 ± 0.17 
0.80 ± 0.29 
1.19 ±0.43 


3.0 (0.08) 
0.1 (0.75) 
1.0 (0.32) 
1.9 (0.17) 
3.0 (0.08) 
2.8 (0.09) 
4.5 (0.03) 
5.0 (0.03) 



dN 
dm 



= < 



if rribreak < m < 100 M G 



(5) 

with r\ = —2.35 (Salpeter slope), a = —3 for the bottom- 
heavy IMFs suggested by van Dokkum & Conroy (2010) or 
a — —3.5 in the most extreme case. 

Our results for the bottom-heavy stellar mass fraction 
are listed in Table |4j We also list the likelihood ratios be- 
tween the nominal isotropic Hernquist model lensing and 
dynamic stellar mass fraction and the stellar mass fractions 
obtained from stellar populations and colors for the different 
IMFs (equivalent to Ax 2 if the distribution was Gaussian, 
which we assume as a first approximation), comparing their 
maximum a poteriori difference with the no-difference hy- 
pothesis. Vertical arrows in Figure [5] show the stellar mass 
fractions predicted by stellar population synthesis models 
and SDSS color by these different IMFs and the stellar 
mass fraction obtained with lensing and dynamics. Using 
w.break = 0.7Mq and a bottom-heavy IMF, we find that 
inferred stellar mass fraction exceeds unity for a = —3.5, 
inconsistent with the lensing mass. An extreme 'bottom- 
heavy' IMF with slope of a = —3.5 is inconsistent at the 
> 90 per cent confidence level with the lensing and kinematic 
results. The a = —3 model is only marginally consistent for 
Tibreak = 0.3 but is also excluded at the > 90 per cent confi- 
dence level for mb roa k = 0.7. A Salpeter IMF gives the best 
agreement with the lensing and kinematics. It is important 
to mention that we do not include any possible effects of 
large-scale structure line-of-sight contamination (e.g. from a 
group elongated along the line of sight), that would decrease 
the total mass assumed here. Recalling Treu et al. (2009) 



and Guimaraes and Laerte Sodre Jr (2011), SLACS lenses 
are shown to be unbiased sample in relation to a random 
LOS, despite the fact that the lenses are elliptical galaxies 
which are often found in dense regions. Moreover, other pos- 
sible explanations, such as a top heavy IMF with more black 
holes and neutron stars remnants, are still possible. As dis- 
cussed in van Dokkum (2008), top-heavy IMFs have fewer 
low-mass stars than a standard Salpeter IMF but many more 
high-mass stars. Nevertheless, it is important to clarify that 
since our method only infers the overall M/L, we cannot 
distinguish a Salpeter IMF from a bottom-light IMF like a 
Chabrier IMF due to the presence of remnants (Treu et al. 
2010 and Auger et al. 2010). 



6 SUMMARY & CONCLUSIONS 

In this paper we present the first results from a new spec- 
troscopic survey of massive early-type lens galaxies: The 
X-Shooter Lens Survey (XLENS). The combination high- 
fidelity UVIS-IR spectroscopy from the X-Shooter instru- 
ment on the VLT, with the strong gravitational lensing mass 
determination has enabled us to conduct an in-depth study 
of the internal structure of the luminous elliptical galaxy 
SDSS J1148+1930 at z = 0.444. We find the following: 

(i) The luminosity-weighted stellar velocity dispersion is 
(a*)(< -Reff) = 352 ± 10 ± 16kms _1 , more accurate and 
considerably lower than a previously published value of ~ 
450 km s -1 . 

(ii) A single-component (stellar plus dark) mass model of 
the lens galaxy yields a logarithmic total-density slope of 
i = 1.72±g;g| (68 per cent CL; p tot oc r~ 7 '). 

(iii) The projected stellar mass fraction, derived solely 
from the lensing and dynamical data, is /*(< -Re) = 
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0.19^0 09 (68 per cent CL) inside the Einstein radius for a 
Hernquist profile and no anisotropy. The dark-matter frac- 
tion inside the effective radius /dm(< Res) = 0.60+q gg±0.1 
(68 per cent CL), where the latter error is systematic. 

(iv) Based on the SDSS colors, we find /»,Sai P (< Re) = 
0.17 ± 0.06 for a Salpeter IMF and /*, C hab(< Re) = 0.07 ± 
0.018 for a Chabrier IMF. A Salpeter IMF gives the best 
agreement between lensing and dynamics constraints on the 
stellar mass fraction, therefore it is preferred to a Chabrier 
IMF. Dwarf-rich IMFs with a > 3 (with dN/dM oc M~ a ) 
in the lower mass range of 0. 1-0.7 M@, - such as those re- 
cently suggested by van Dokkum & Conroy (2010) for mas- 
sive early-type galaxies (a = —3) - are excluded at the > 90 
per cent CL and in some cases (a = —3.5) violate the total 
lensing-derived mass limit. 

This massive early-type galaxy lies at the extreme end of the 
trend found by, e.g., Auger et al. (2010) that the dark matter 
fraction within the effective radii is a monotonically increas- 
ing function of galaxy mass. In fact SDSS J1148+1930 is 
already dark-matter dominated within that region. We find 
that a Salpeter IMF agrees best with the total stellar mass 
derived from lensing and stellar kinematics as well as with 
its SDSS colors. As in Treu et al. (2010) and Grillo & Gobat 
(2010), a Salpeter IMF appears to be the best option for very 
massive early-type galaxies. Although slightly more massive 
IMFs cannot be excluded given the typical uncertainties. 
A bottom-light IMF such as a Chabrier (or Kroupa) IMF 
agrees only marginally and we exclude a steep 'dwarf-rich' 
IMF with a = —3.5 at > 90 per cent CL. Somewhat shal- 
lower IMFs with a ~ —3.0, as suggested by van Dokkum 
& Conroy (2010), are marginally acceptable. We conclude 
that our data are fully consistent with SDSS J1148+1930 
being a massive early-type galaxy that is dark-matter dom- 
inated inside its effective radius and having a Salpeter IMF. 
No strong evidence for an even more bottom-heavy IMF is 
found, consistent with previous results (Treu et al. 2010), 
although uncertainties are still large.. 

Further studies are required to break the degeneracy 
between the central dark-matter fraction and distribution 
and the stellar IMF. In forthcoming papers of the XLENS 
survey, we will extend the study to more massive systems at 
z > 0.5 and also a sub-sample of SLACS lenses at z < 0.5. 
Observations are ongoing. In those papers, we will perform 
more detailed stellar population analyses using the full UV- 
optical-NIR spectrum and obtain independent contraints on 
the the physical parameters that may correlate with IMF 
normalization (i.e., age and metallicity) or that may be the 
cause of the correlation between dark-matter content and 
velocity dispersion. 
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